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We search for viable f(R) theories of gravity, making use of the equivalence between such theories and 
scalar-tensor gravity. We find that models can be made consistent with solar system constraints either by giving 
the scalar a high mass or by exploiting the so-called chameleon effect. However, in both cases, it appears likely 
that any late-time cosmic acceleration will be observationally indistinguishable from acceleration caused by a 
cosmological constant. We also explore further observational constraints from, e.g., big bang nucleosynthesis 
and inflation. 



I. INTRODUCTION 

Although the emerging cosmological standard model fits 
measurements spectacularly well (see d 01 for recent re- 
views), it raises three pressing questions: what is the physics 
of the postulated dark matter, dark energy and inflation en- 
ergy? The need to postulate the existence of as many as three 
new substances to fit the data has caused unease among some 
cosmologists JH 0, HI |a] and prompted concern that these 
complicated dark matter flavors constitute a modern form of 
epicycles. Our only knowledge about these purported sub- 
stances comes from their gravitational effects. There have 
therefore been numerous suggestions that the apparent com- 
plications can be eliminated by modifying the laws of grav- 
ity to remove the need for dark matter J3, @], dark energy 
ia OH [Till and inflation lfl2ll . and perhaps even all three to- 
gether [13]. Since attempts to explain away dark matter with 
modified gravity have been severely challenged by recent ob- 
servations, notably of the so-called bullet cluster 11411 . we will 
focus on dark energy (hereafter "DE") and inflation. 

There is also a second motivation for exploring alternative 
gravity theories: observational constraints on parametrized 
departures from general relativity (GR) have provided increas- 
ingly precise tests of GR and elevated confidence in its validity 

diii. 



A. f(R) gravity 

An extensively studied generalization of general relativ- 
ity involves modifying the Einstein-Hilbert Lagrangian in the 
simplest possible way, replacing R — 2 A by a more general 

function /(E) 1 d M, IM Hi HI HI HI m m lli m 



For the general case where / depends on the full Riemann tensor R^ a a 
rather than merely on its contraction into the Ricci scalar R, this program 
is more complicated; a subset of these theories which are ghost free can 
be written as f(R, G), where G = R \ llJa R g ^ l ' a - iR^ R^ + R 2 
is the Gauss-Bonnet scalar in 4 dimensions 11711 . These theories lack a 
simple description in terms of canonical fields; there is no so-called Ein- 
stein Frame. Progress has nevertheless been made along these lines, and 
such Lagrangians may have more relevance to DE fl7l [lil fl^l than ones 
independent of G. 



l30l I3TI l32tl . The equations of motion derived from this La- 
grangian differ from Einstein's field equations when f(R) is 
nonlinear, but the theory retains the elegant property of gen- 
eral coordinate invariance. In such a theory, the acceleration 
of our universe may be explained if f(R) departs from lin- 
earity at small R, corresponding to late times in cosmological 
evolution. In this case it may be possible to avoid invoking 
a cosmological constant to explain cosmic acceleration, al- 
though one then replaces the problem of a small cosmological 
constant with the problem of no cosmological constant. In 
such models, the effective DE is dynamic, i.e., it is not equiv- 
alent to a cosmological constant, leading to potentially inter- 
esting observational signatures. We refer to these models as 
f(R) -DE theories. 

In addition to potentially explaining late-time acceleration, 
f(R) theories may be relevant to early-universe physics, par- 
ticularly if f(R) is non-linear at large R 1 12113311 . More gener- 
ally, it is of interest to study f(R) theories because they are ar- 
guable the simplest setting in which one can attack the general 
question of which modified theories of gravity are allowed. 



B. The equivalence with scalar tensor gravity 

The modified Einstein field equations (and so the new 
Friedmann equation) resulting from a non-linear f(R) in the 
action can be seen simply as the addition of a new scalar de- 
gree of freedom. In particular, it is well-known that these the- 
ories are exactly equivalent to a scalar-tensor theory |[34l[3~5ll . 
It is therefore no surprise that for f(R)-DE theories, it is this 
scalar which drives the DE. Before reviewing the mathematics 
of this equivalence in full detail in section [II] we will discuss 
some important qualitative features below. 

One can discuss the theory in terms of the original met- 
ric g^v, in which case the degrees of freedom are not man- 
ifest. Alternatively, by a conformal relabeling, one can re- 
veal the theory to be regular gravity g^ v , plus scalar field <\>. 
The former viewpoint is referred to as the Jordan Frame (JF) 
and the latter as the Einstein Frame (EF). Here qb has the pe- 
culiar feature that in the JF, it exactly determines the Ricci 
scalar R and vice versa. So in the JF, the Ricci scalar can 
in a sense be considered a non-canonical yet dynamical scalar 
field. This feature is absent in normal general relativity, where 
R = -T/M^y + 4A is algebraically fixed by the trace of the 
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stress energy tensor T. Working in either frame is satisfactory 
as long as one is careful about what quantities are actually 
measurable, but we will find that the EF is much more useful 
for most of our calculations. 

The coupling of the scalar field to matter is fixed in f(R) 
gravity, and is essentially of the same strength as the cou- 
pling of the graviton to matter, except for the important case 
of massless conformally invariant fields, which do not couple 
to <j> at all. The dynamics of the theory are completely speci- 
fied by the potential V(<f>) for the scalar field in the EF, which 
is uniquely determined by the functional form of f(R). 

C. The R - n 2 /R example 

Such a scalar field is not without observational consequence 
for solar system tests of gravity, especially for /(i?)-DE mod- 
els. For any scalar field driving DE, we can come to the fol- 
lowing conclusions: First, the field value must vary on a 
time scale of order the Hubble time Hq , if the DE is distin- 
guishable from a cosmological constant (for longer time scale, 
the DE looks like a cosmological constant; for shorter time 
scales, we no longer get acceleration). On general grounds, 
such a scalar field must have a mass of order rr& ~ Hq. Sec- 
ond, the Compton wavelength of this scalar field is on the or- 
der of the Hubble distance, so it will mediate an attractive fifth 
force which is distinguishable from gravity by the absence of 
any coupling to light. Unless the coupling to matter is tiny 
compared to that of gravity, many solar system based tests of 
gravity would fail, such as measurements of the bending of 
light around the Sun lfl5i[l6ll . 

The archetypal example of f(R)-DE suffers from problems 
such as these. This model invokes the function [ l36ll . 

f(R) = R-£ (1) 

for /it w Hq. This gives a V((f>) in the EF with 
a runaway exponential potential at late times: V(<f>) ~ 
^o-^pi ex P( — (3/V6)^/Afpi) (here large <f> means small R 
which means late times.) With no matter, such a potential in 
the JF gives rise to an accelerating universe with the equation 
of state parameter % = —2/3 03611 . This model, however, 
is riddled with problems. First, the theory does not pass so- 
lar system tests l20l l2ll I22I l30ll . and second, the cosmology 
is inconsistent with observation when non-relativistic matter 
is present l23ll . Both problems can be understood in the dual 
scalar tensor theory. 

For cosmology, during the matter dominated phase but at 
high redshifts, the influence on the dynamics of <f> from the 
potential V is small compared to the influence from the cou- 
pling to matter, which manifests itself in terms of an effective 
potential for </> of the form 

y cff (0) = v(4>) + p NR exp . ( 2 ) 

where Pnr is the energy density of non-relativistic (NR) mat- 
ter. (More details of the exact form of this potential will be 



presented in the next section.) The second term dominates 
because HqM^ <C pnr and (f> then rolls down the poten- 
tial generated by pnr and not V. The result is that the uni- 
verse is driven away from the expected matter dominated era 
(MDE) into a radiation dominated expansion in the JF with 
H 2 oc a -4 , after which it crosses directly into the acceler- 
ating phase, with expansion driven by DE with an effective 
equation of state parameter 11; = —2/3. This special radiation- 
dominated-like phase (which is not driven by radiation) was 
dubbed the </>MDE by 112 311 . where it was made clear that this 
phase is inconsistent with observation. We say that this poten- 
tial V is unstable to large cosmological non-relativistic densi- 
ties. 

For the solar system tests, the potential V(<p) is also negli- 
gible, so the theory behaves exactly like a scalar-tensor theory 
with no potential. Because the coupling to matter has the same 
strength as that to gravity, the scalar field mediates a long- 
range fifth force, and the theory is ruled out by solar system 
tests. In particular, i22ll found that 7 = 1/2 in the PPN frame- 
work, which is in gross violation of the experimental bound. 

The above solar system tests also seem to rule out more 
general classes of /(i?) -DE models (TUMHSIHIll. How- 
ever on the cosmology front, it seems that one can cook up 
examples of f(R) consistent with some dynamical dark en- 
ergy 1241 125, 28]: by demanding that the cosmological expan- 
sion a(t) take a certain form, one can integrate a differential 
equation for the function / that by design gives a universe 
with any desired expansion history a(t). In this way, one gets 
around the cosmological instability of the archetypal model 
mentioned above. However, these functions are arguably very 
contrived, and further investigation of solar system predictions 
are required to determine whether these models are viable. 



D. What /(i?) -theories are allowed? 

We now try to find viable f(R) theories by examining what 
is acceptable on the scalar tensor side. We focus on theories 
that pass solar system tests. Because the coupling of the scalar 
field to matter is fixed in f(R) theories, and the only freedom 
we have is with the potential V, we must choose V in such 
a way as to hide the scalar field from the solar system tests 
that caused problems for the models described above. We are 
aware of only two ways to do this. The first is the Chameleon 
scalar field, which uses non-linear effects from a very specific 
singular form of potential to hide the scalar field from current 
tests of gravity 11371 l38tl . The second is simply to give the 
scalar field a quadratic potential with mass > 10~ 3 eV, so 
that the fifth force has an extent less than 0.2mm and so cannot 
be currently measured by laboratory searches for a fifth force 
& 

We will find simple f(R) models which reproduce these 
two types of potentials and so by design pass solar system 
tests. Finding functions / which give exactly these potentials 
will simply generate models which are indistinguishable from 
their scalar-tensor equivalent. However if we search for sim- 
ple choices of / that reproduce these potentials in a certain 
limit, then these theories will not be exactly equivalent and 
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might have distinguishable features. 

The Chameleon type f(R) model seems to be the most 
plausible model for attacking DE, as at first glance it seems 
to get around the general problems mentioned above. Indeed, 
one Chameleon model will arise quite naturally from a simple 
choice of /. However, we will show that the solar system con- 
straints on this model preclude any possible interesting late- 
time cosmological behavior: the acceleration is observation- 
ally indistinguishable from a cosmological constant. In par- 
ticular, for all the relevant physical situations this Chameleon 
model is the same as has been considered before with no dis- 
tinguishing features. However, this model might provide clues 
in a search for viable f(R) theories that pass solar system tests 
and that may give interesting late-time behavior. 

In an independent recent analysis, 13111 also discussed the 
Chameleon effect in f(R) theories. They focus on a slightly 
different set of Chameleon potentials and come to similar 
conclusions. Their results and ours together suggest that the 
Chameleon effect may be generic to f(R) theories. 

We now turn from attempts to explain DE in f(R) mod- 
els to an arguably more plausible scenario, which is simply 
to give the scalar field a large mass. These models have no 
relevance for dynamic DE, but they do have interesting con- 
sequences for early universe cosmology. The most theoreti- 
cally best motivated functions, namely polynomials in R, fit 
this class of f(R) theories. The aim of this investigation is 
to explore what we can possibly know about the function /. 
Because this question is very general, we will restrict our at- 
tention to a sub-class of plausible f(R) models. 

For these polynomial models, we will investigate possible 
inflationary scenarios where the scalar partner <f> is the infla- 
ton. We find the relevant model parameters which seed the 
fluctuations of the CMB in accordance with experiment. We 
then investigate general constraints on the model parameters 
where <j> is n °t an inflaton. We use solar system tests, nu- 
cleosynthesis constraints and finally an instability which is 
present in these theories when another slow roll inflaton ijj is 
invoked to explain CMB fluctuations. This instability is anal- 
ogous to that of the 0MDE described above. 

The rest of this paper is organized as follows. In section 
[II] we review the equivalence of f(R) theories with scalar 
tensor theories, elucidating all the essential points we will 
need to proceed. Then in sections [Til] and [IV] we explore the 
Chameleon model and massive theories, respectively, focus- 
ing on observational constraints. We summarize our conclu- 
sions in sectionlVl 



ceed. See for example fill l35i l40ll for more details of the 
equivalence with scalar tensor theories. 

We choose to fix the connection in R as the Christof- 
fel symbols and not an independent field, as opposed to the 
Palatini formalism - which results in a very different theory 

HUE BIB El • 

If one simply varies the action Eq. (O with respect to the 
metric g M „, then a fourth order equation for the metric results. 
One can argue (using general coordinate invariance) that the 
degrees of freedom in the field g M „ can be split into a massless 
spin-2 field g„ y and a massive scalar field cf> with second order 
equations of motion. This split is easily revealed at the level 
of the action. Following for example 14011 we introduce a new 
auxiliary scalar field Q (a Lagrange multiplier) the gravity part 
of Eq. (01 may be written as 



5grav — 



-g-^ (f(Q)(R- 



Q) + f(Q)) (4) 



As long as f"(Q) 0, the equation of motion (S/SQ) gives 
Q = R and Eq. becomes the original gravity action. This 
may be written in the more suggestive form 




x 2 v( x ) 



(5) 



by relabeling f'(Q) = X- This is a scalar tensor theory of 
gravity with Brans Dicke parameter lobd = [46] and poten- 
tial 



v(x) = ^[Q(x)x-f(Q(x))} 



(6) 



Here Q(x) solves \ — f'(Q)- Finally a rescaling of the metric 
(which should be thought of as a field relabeling) 



reveals the kinetic terms for the scalar field: 



(7) 



Sep = J d**V=I (^R - ^d^d^ - 



(8) 



where the new canonical scalar field <fi is related to x, Q, R 
through 



n. f(R) DUALITY WITH SCALAR TENSOR THEORIES 

We study the "modified" gravity theory defined by the ac- 
tion 

S JF = J d 4 x^5^/(i?) + S M [9^ *, A a , . . .] (3) 

Here we present a run down of the map to the scalar tensor 
theory; displaying the most important points needed to pro- 



f'(R) = f'(Q) = X 



exp i 



(9) 



As the kinetic terms for g^ v and <f> are now both canonical, we 
see that these are the true degrees of freedom of f(R) grav- 
ity. This demonstrates that the theories defined by Sjf ( the 
Jordan Frame) and Sef (the Einstein Frame) are completely 
equivalent when f"(Q) ^ 0. We choose to analyze the theory 
in the Einstein Frame as things are much simpler here. It is, 
however, important to be careful to interpret results correctly, 
making reference to what is observed. In particular, matter is 
defined in the Jordan Frame, and hence it will be most sensible 
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to talk about JF observables. We will give a simple example 
of this when we have introduced some matter. 

The equations of motion for <j> resulting from Eq. I© are 



V6M pl ' 



(10) 



and for the metric g^, 

R,u - \g„ v R = A/p? (t# + f *,) (11) 

with the energy momentum tensors 

T$ = X^Tll U' 1 ^ ■ ■ •) (12) 
f*, = dpWut + 9^ [~\9 ap d a( t>d^ + V(cf>^ (13) 

Note that only the combination T~ + is conserved in the 
EF. 

There are two important observations to be made about 
Eq. (O relating to the extra coupling to matter. First, the 
T M /Mpis/E term in Eq. ( TTOb represents an additional density 
dependent "force" on the scalar field, and for special cases 
where we can solve for the qb dependence of T M /M p \y6 ex- 
plicitly, as in 13811 . we can think of the scalar field living in an 
effective potential. We will see two examples where this force 
is important, the most dramatic being the Chameleon effect. 

Second, <\> couples to matter as strongly as conventional 
gravity {g^) does. Hence, as was already mentioned, <j> will 
mediate a detectable fifth force for solar system tests unless 
we do something dramatic to hide it. Finding theories which 
hide <f> from solar system tests is the focus of this paper. 



A. Matter and Cosmology in f(R) theories 

Let us first consider the coupling to standard model fields, 
assuming that they are defined in the JF. This is important for 
understanding how (j> may decay. Massless scalar fields con- 
formally coupled to gravity and massless gauge bosons be- 
have the same in the two frames and so do not couple to <fi. 
However, a minimally coupled (real) scalar field <£> and a Dirac 
field if? have extra interactions with 6 in the EF: 



S'<i 



d xy —g 



■ l(d$Y -^mlx- 1 ® 2 (14) 



1 



V6M, 



-<f>g^d^d u , 



pi 

1/2 W. 



(15) 



where the JF fields have been rescaled as $ = x _1 / 2 $ and 
if? = x~ 3 ^ 4 ^- Note that the cosmologically evolving field 
<j) = 4>(t) will change the masses of the standard model parti- 
cles in the EF as 



m 



mx 



-1/2 



(16) 



and small excitations 5(j) around the average value 4>{t) will 
roughly speaking interact via the vertices defined by the inter- 
action Lagrangian, 



1 



™2 



\/6M 



m^S^f) (17) 



pi 



to lowest order in 1/M P \. The mass shift in Eq. (TToT l has an 
interesting consequence in the EF; it shifts the frequency of 
the absorption and emission lines by a factor of x -1 ^ 2 - This 
effect will be indistinguishable from the normal cosmological 
redshift due to expansion, and our effective redshift will be the 
combination of both cosmological expansion and mass shift: 
(1 + z)~ x = dx -1 / 2 , where a is the scale factor in the EF 
normalized equal to unity today. This combination turns out to 
be the Jordan frame scale factor a (see below), so our redshift 
measurements coincide in both frames as expected. 

Perfect fluids are best examined in the JF, because it is here 
that their energy momentum tensor is conserved. For a gen- 
eral JF metric one can solve for the flow of the fluid using 
conservation of T^l and number flux nlf^ (or other relevant 
physical principles) and then map into the EF via Eq. ( fT2l . 



1. The homogeneous and isotropic case 

For example, consider a homogeneous isotropic cosmol- 
ogy, 

(JF) ds 2 = dt 2 - a(t) 2 di 2 , U" = {dtY , (18) 
(EF) ds 2 = dt 2 - a(t) 2 di 2 , U" = {dif , (19) 

where C/ M and C/ M are the local fluid velocities in the two 
frames. The quantities above are related by 

d = X ^a, dt = x^dt, U» = X ~ 1/2 U^ . (20) 
These relations imply that 



H = y^ 2 H- ^ 



V6M, 



(21) 



pi 



For example, applying the principles of entropy "conserva- 
tion" and number conservation in the JF (one may also need 
to demand thermal and chemical equilibrium as relevant to the 
early universe) results in known p(a) and p(a) such that the 
EF energy momentum tensor may be written as 



fj* = pU^U v + p { U^U V + g^ 



where 



P = X P 



(V 1/2 ) 



-2 /- -1/2 

P = X P(a>X 



(22) 



(23) 



The cosmological equations of motion are, 



3H 2 M 2 =p+^(a t 0) +V{<f>), (24) 
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d 2 cj)+3Hd t d = - 



dV eR {4>,a) _ dV E T 



M 



06 



V6M, 



(25) 



pi 



The effective potential for the scalar field coupled to homoge- 
neous and isotropic matter is 

V oS {<t>, ~a) = V{6) + p = V{4>) + X ~ 2 P (a X - 1/2 ) . (26) 

For the special case where the only density components 
present are non-relativistic (p = pnr oc a~ 3 ) and ultra- 
relativistic (p = pr oc a~ 4 ) fluids, the effective potential is 



V e g(d) = V(6)+pMa)e M * vs +PR(a) (27) 

where for convenience we define Pnr(o) = 
X _3/2 PNR(ax _1/2 ) oc &- 3 and p K (a) = x~ 2 PR(ax~ 1/2 ) oc 
a -4 . Note that relativistic particles provide no force on <f> 
because T vanishes, or equivalently because Pr(o) appears 
simply as an additive constant to the potential in Eq. ([27). 



2. The spherically symmetric case 

We now turn to the case of a spherically symmetric distri- 
bution of non-relativistic matter pnr(^) in the JF, for which 
we aim to solve for the metric <? M „. We wish to consider this 
problem in the EF, where & will take a spherically symmetric 
form and gravity behaves like GR coupled to p = x~ 2 Pnr- 
In the weak field limit, we write the metrics in the two frames 



(JF) ds 2 = -{l-2A(r))dt 2 + (l + 2B(r))dr 2 

+ r 2 dfl 2 (28a) 

(EF) ds 2 = -{1 - 2A(r))dt 2 + (1 + 2B(r))df 2 

+ r 2 dQ 2 (28b) 

where r — x r and for small <j)/M p i, the gravitational po- 
tentials are related by 



A(r) A(f) + 
B(r) « B(r) + 



1 #(f) 
V6M p i dlnf 



(29a) 
(29b) 



Following 13811 we define a non-relativistic energy density 
Pnr(^) = X~ 3 ^ 2 p( r ) m me EF which is conserved there and 
is analogous to Pnr(S) defined above for cosmology. Ignor- 
ing the back reaction of the metric on 6, we take i=s -q^ 
in Eq. (|T0l> and find as in [38] that 



1 d ( _ 2 dd 



_iPnr(/0 dVeg(<l>,r) 



r 2 dr V drj V ' {<P) X ^QM pl 
where again the effective potential is 

Kff = V{4>) + x~ 1/2 pNR(r) ■ 



(30) 
(31) 



Solving Eq. d30b for 6 then allows us to find the metric in the 
JF via Eq. ([29). 

As an instructive example, consider the quadratic potential 
V(<p) — ni^d 2 /2 and a uniform sphere of mass M c and radius 
R c , the solution external to the sphere is given by a Yukawa 
potential 



1 M r e 



M, 



p i V6 47rM 2 !r 



(32) 



assuming that m^R <C 1 and (j)/M p \ <C 1 so that f w r. If 
we ignore the energy density of the profile 4>(r), then outside 
the object there is vacuum. The metric in the EF is then sim- 
ply the Schwarzschild solution for mass M c . In other words, 
the potentials in Eq. ( I28at are given by A(r) = B(r) — 
M c /8wM^r in the weak field limit \A\, \B\ < 1. In the 
JF using Eq. (f29), one finds the corresponding potentials 



A(r) w A(r) 
B(r) w A(r) 



1+ r 



1 - -e" 
3 



"(1 + m^r) 



(33a) 



(33b) 



For r 1 we find that the PPN parameter 7 = 1/2, a well 

known result for a Brans Dicke theory j46Tl with ujbd = 
0. 

The key feature here is the effective potential from Eqs. ( f27b 
and PP . We have now seen that it makes a crucial difference 
in two situations, and it will play an important role in the next 
two sections as well. 



in. AN f(R) CHAMELEON 

In this section, we consider f(R) theories that are able to 
pass solar system tests of gravity because of the so-called 
"Chameleon" effect. We first present a theory that is by de- 
sign very similar to the original Chameleon model presented 
in iTsvIl . We will give a brief description of how this model 
evades solar system constraints, and then move on to the cos- 
mology of these f(R) theories, concentrating in particular on 
their relation to DE. Throu ghou t this discussion we refer the 
reader to the original work B37ll38ll47ll . highlighting the dif- 
ferences between the original and f(R) Chameleons. 

The Chameleon model belongs to the following general 
class of models, 



f(R) = R-(1- m)^ 2 



2 A. 



(34) 



The sign of the second factor is important to reproduce the 
Chameleon, and the (1 — m) factor ensures that the theory is 
equivalent to GR asnn 1. These models have been consid- 
ered before in the literature fill, l23ll : in particular, this class 
contains the original DE f(R) of Eq. ([T) when m = — 1 and 
A = 0. 

The potential for 6 in the EF is 



V&) 



2x 2 



(m - l) 2 



x-i 



M 2 ,A 



(35) 
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where \ = exp(i/2/3(/>/Mpi) as usual, 
for |0/Mpi| <C 1, this reduces to 



For < m < 1 and 



M p 2 ,A, 



(36) 



defined for <f> < 0, where the old parameters //, m are related 
to the new parameters M, n through 



1 



M 2 = 



(2(1 



M 4+r 



n+2 



(37) 



pi 



The preferred values used in ll3~7ll are M ~ 10~ 3 eV and n ~ 
1. In the f(R) theory, these values give m ~ 1/2 and /i ~ 
10 _50 eV, i.e. much smaller than the Hubble scale today. 

For small |0|/M p i, this singular potential is equivalent to 
the potential considered in [ 37 ] for the Chameleon scalar field, 
albeit with (f> — > — (f>. The coupling to matter, which is a very 
important feature of this model, is also very similar. In 13711 . a 
species of particles i is assumed to have its own Jordan Frame 
metric gfy, with respect to which it is defined and a conformal 
coupling to the metric in the EF 



(i) _ 2/3,<p/M pl ~ 

yuv ~ c 



(38) 



Comparing this to Eq. (0, the f(R) Chameleon has /3, = 
— l/v6 for all matter species, so that all the Jordan Frame 
metrics coincide. 

In the original Chameleon model, the fli were specifically 
chosen to be different so that (j> would show up in tests of 
the weak equivalence principle (WEP). The f(R) Chameleon 
does not show up in tests of the WEP, so the solar system 
constraints will be less stringent here. 

This coupling to matter, along with the singular poten- 
tial, are the defining features of this f(R) that make it a 
Chameleon theory. The effective potential V c s, discussed in 
the previous section (see for example Eq. |27| ). is then a bal- 
ance between two forces; V pushing <f> toward more negative 
values and the density-dependent term pushing <j> toward more 
positive values. So although the singular potential Eq. d36t 
has no minimum and hence no stable "vacuum", the effec- 
tive potential Eq. d2Tb including the coupling to matter does 
have a minimum. In fact, the density dependent term pushes 
the scalar field <f> up against the potential wall created by the 
singularity in V at <f> = 0. Indeed, the field value </> m i n at 
the minimum of the effective potential V c s and the mass 
of 0's excitation around that given minimum are both highly 
sensitive increasing functions of the background density pnr, 
as illustrated in Fig.Q] Using Eq. OTb for small \<j>\/M v \, the 
field value at the minimum and the curvature of the minimum 
are, respectively, 



m(l - m) ( 



\/6M p i 



2 
2 



Pnr 



Pnr 



V6M 



pi 



3(1 - m) M 



(39) 
(40) 



pi 



It is plausible that a scalar field <p which is very light for cos- 
mological densities is heavy for solar system densities and 
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FIG. 1: Effective potential for the Chameleon model Eq. d35 b with 
decreasing pnr/ fJ 2 = 100, 50, 20 and m = 1. Note that (ftmin 
and the mass (the curvature of the minimum) are very sensitive 
to the background energy density pnr. 



hence currently undetectable. However, as we will now see, 
the actual mechanism that "hides" <f> from solar system tests is 
a bit more complicated than this. 



A. Solar System Tests 

In this section, we will derive solar system and laboratory 
constraints on the parameters (fi,m), summarized in Figure 
|2] The profile of <f>(f) in the solar system (around the Earth, 
around the Sun, etc.) is of interest for solar system tests of 
gravity: it determines the size of the fifth force and the post- 
Newtonian parameter 7. Because the effective potential for 
<j> changes in different density environments, the differential 
equation governing the profile (j>(f ) in Eq. (f30b is highly non- 
linear, and the standard Yukawa profile of Eq. d32b does not 
always arise. These non-linear features have been studied in 
|38y, where it was found that for a spherically symmetric ob- 
ject of mass M c and radius R c surrounded by a gas of asymp- 
totic density p^, the profile is governed by the so-called "thin- 
shell" parameter, 



\4>l 



2AttM^R c 



pi 



Mr. 



(41) 



where 0^ in and ^>^ lin are the minima of the effective potential 
in the presence of the asymptotic energy density = 
and central energy densities /?nr = Pc respectively, see 
Eq. (1391 . If A is large, then the external profile of cj> is the 
usual Yukawa profile Eq. ( l32t with mass = (pnr = 
Poo), the curvature of the effective potential in the presence of 
the asymptotic density; see Eq. (1401 . If A is small, then the 
Yukawa profile is suppressed by a factor of A. The term "thin 
shell" comes from the fact that only a portion of such a "thin 
shell" object contributes to the external Yukawa profile, the 
thickness of the shell being roughly (Ai? c ). We simply treat 
A as a parameter that suppresses this profile if A <C 1. 



7 



For example, let us consider the profile <f> around the Sun, 
with M c = Mg un and R c = i?sun- Assuming that we are in 
the thin shell regime (A <C 1), the Yukawa profile of Eq. < 
suppressed by a factor A becomes, 



-) 



A M Sun e- 



(42) 



V6 4^Af pl r 

As in [37], we take the asymptotic density used to find <poo 
and moo as that of the local homogeneous density of dark and 
baryonic matter in our Galaxy: p^ s» 10 _24 g/cm 3 . Follow- 
ing the discussion in Section [II] the metric in the EF external 
to the Sun is just the Schwarzschild metric (in the weak field 
limit) with Newtonian potential A(r) w Ms un / '(8-irM^r). 



Using Eq. (|29j l to map this metric into the JF metric g^ v 
X _1 <W' we find 



ds 2 



I - 2 Sir i ( 1 + je 



A 



r 2 dn 2 



1 + 2A(r) ( 1 - — e~ m ~ r (l + m^r) 



dr\AJ) 



Assuming that the Compton wavelength m^ 1 is much larger 
than solar system scales (we will confirm this later), we obtain 
within the PPN formalism [15] that 
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3- A 
3 + A 



1- (2/3)A 



(44) 



There are several observational constraints on \y — 1|, in- 
cluding ones from the deflection of light and from Shapiro 
time delay. The tightest solar system constraint comes from 
Cassini tracking, giving I7 — 1| < 2.3 x 10~ 5 [16]. Thus 
the "thin shell" parameter satisfies A < 3.5 x 10~ 5 . We take 
l^minl ^ I '/'mini because of the sensitive dependence on the 
local density, so the definition of A in Eq. ( f4Tb becomes 



A « 3|0£n|/V6M p ii(fl 



Sun J 



(45) 



where A(R C = i?s U n) ~ 10~ 6 is the Newtonian potential at 
the surface of the Sun. Using Eq. ( f39l > with /5nr = Poo ~ 
10~ 24 g/cm 3 gives the constraint 



Hi 



<3 



m(l — m) 



-6-Hn 



(46) 



on the theory parameters p and m. For theories which fail 
this bound, we find that the Compton wavelength of tfi for 
the asymptotic background density of our galaxy satisfies 
m oo ^ 10 10 AU. This confirms the assumption that m^ 1 
is large compared to solar system scales, which was used to 
derive this bound. 

As was already noted, the solar system constraints derived 
in 138] are more restrictive. This is because they demanded 
that the couplings {j3i) to different species of particles in equa- 
tion ( f38l > be different. This gives violations of the weak equiv- 
alence principle on Earth-based experiments unless the Earth 
and atmosphere have a thin shell. However, in the f(R) 
Chameleon model, all the /3j are the same, so there will be 
no weak equivalence principle violations. 



The f(R) Chameleon may still show up in searches for a 
fifth force, in particular in tests of the inverse square law. The 
strongest comes from Earth-based laboratory tests of gravity 
such as in the Eot-Wash experiments [39]. By demanding that 
the test masses acquire thin shells, [ 38] found constraints on 
the parameters (M, n) which map into the following bound 
on the /(i?)-parameters (p, m): 



H 2 ~ 11 



m) 



m(l — m) 



1(T 



(47) 



B. Cosmology 

We now turn to the cosmology of the Chameleon scalar 
field, which was studied in 14711 . It was found there and al- 
ready commented on in l38ll that the mass of <fi on cosmolog- 
ical scales is not small enough to give any interesting DE be- 
havior for M 10~ 3 cV and n ~ 1. We will revisit this ques- 
tion in the f(R) context: do any allowed parameters (p,m) 
in Eq. (l47T i give non-vanilla DE? Will there be any cosmolog- 
ically observable differences between this f(R) Chameleon 
and the original model (which is in principle possible because 
higher order terms in the expansion of V in Eq. d3~5T > may 
become important)? We will see that the answer to both of 
these questions is no for the same reason: solar system tests 
preclude the minimum of the effective potential from lying at 
(f> < — M p i on cosmological scales today. 

Let us try to understand this by looking at the details of 
Chameleon cosmology. We first note that, as opposed to l47tl . 
we do not fix AA/^ = M 4 , so we are less restricted as to 
what M can be. The essence of the argument, however, is the 
same as in [47]. Working in the EF, for a large set of initial 
conditions in the early universe, <f> is attracted to the mini- 
mum of the effective potential given by Eq. d27l . The scalar 
field tracks the minimum, which shifts <f)(d) = m i n as the 
universe expands. The energy density in coherent oscillations 
around this minimum are negligible and so there is no "mod- 
uli problem". ( In contrast, this may be a problem for the case 
considered below in Section llVl ) 

We will see that the condition for such a tracking solution 
to be valid is that the minimum satisfies 



(a)/M pl < 1, 



(48) 



so we consistently make this assumption to derive properties 
of the tracking minimum. After matter-radiation equality we 
have the tracking solution 



(a) m(l — m) 



1 — m 



V6M, 



pi 



/0 NR (a) + 4V((j)(a)) 



(49) 

Along this tracking solution, the curvature (mass) around the 
minimum and the speed of the minimum are, respectively 



V6M pl \ fp NR {a)+iV^(d)Y 



H 2 ~ 1 - m \ -0(5) J V Pnr(5) + V(^(o)) J ' 

(50) 
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-1 d<f>{a) 
M pl H dt 



= -3 



-4>(a)\ (l-m)p NR (a) 



M P i J PNR(a) + 4V(^(5)) 



(51) 



Since <fi will track the minimum while m^(5) ^> -ff, Eq. (f50b 
shows that the assumption of Eq. (l48l l is indeed consistent. 

During radiation domination, m^(a)/H 2 ~ 
(— M p i/0(5)) o/omr, so it is possible that at early times the 
scalar field is unbound at early times. We know the expansion 
history and the effective value of Newtons constant Gn quite 
well JMH around BBN; if <p is unbound, we have no reason 
to believe that Gat is near today's value. Requiring that it is 
bound before the beginning of BBN gives a constraint that we 
have included in Figure [2] 

Returning to the matter-dominated era, Eq. (l48T > implies that 
the expansion history in the JF may be written as 



Pnr(o) + V(0(ao)) 



(52) 



For \<fi(ao)\/M p \ <C 1 today, this is just the usual Fried- 
mann equation with a cosmological constant, where in ac- 
cordance with experiment we are forced to identify V(<p(ao)) 
with px (0), the current dark energy density. Note that the pa- 
rameter A in V, which we have not fixed, allows us to make 
this choice independent of any values of p and to. For m not 
small, AMpj « px (0); however, for small m we will see later 
that the situation will be slightly different. 

This implies that the only way to get interesting late-time 
cosmological behavior is to not have \(p(ao)/M p \\ <C 1 but 
rather \<f>(ao) /M p \\ ~ 1 today. In this case the tracking so- 
lution above is not valid; the scalar field is no longer stuck 
at the minimum, and we might not have to invoke a constant 
A in V to explain todays accelerated expansion. Rather the 
acceleration would be driven by a quintessence type phase. 

However, one can show that given the solar system con- 
straints, \(f>(ao) / M p \\ ~ 1 is not possible. In fact, as we will 
now show, a stronger statement can be made: even if we con- 
tinue to assume Eq. d48l l. so that the tracker solution is still 
valid, the solutions that are consistent with solar system tests 
always give DE behavior that is "vanilla," i.e., indistinguish- 
able from a cosmological constant. 

In these models, the effective dark energy density is 



Px(a) 



V(4>(a)) + (^^j (PNR(a) +px(0)) 



6p N R(fl)(l -m) . 
PNR(a)+4 Px (0)y k } 

where V{<j)(a)) - p x {0) = 0{(j>/M v y). If we expected 
Eq. ( T53l l to give interesting behavior in the allowed region of 
parameter space, we would fit the Friedmann equation with 
px (a) to the combined knowledge of the expansion history 
and find the allowed values of (/i, to). We will instead adopt 
a simpler approach, defining "non-vanilla DE" through the ef- 
fective equation of state parameter, 



w x 



1 din px{o) 
3 din a 



1. 



(54) 



This is the relevant equation of state that one would measure 
from the expansion history (that is not p^/p^,). We say that the 
DE is non-vanilla if \wx + 1| > .01, which is quite-optimistic 
as to future observational capabilities B49I1 . However, because 
our result is null the exact criterion is not important. 

The resulting constraint on p and m is shown in Figure [2] 
along with the solar system constraints. As the Figure shows, 
all models consistent with solar system tests are "vanilla" - 
that is, indistinguishable from a cosmological constant. 




FIG. 2: Solar system constraints on the f(R) Chameleon are seen to 
exclude all models where the "dark energy" is observationally dis- 
tinguishable from a cosmological constant (labeled "dynamic DE"). 
The two different solar system constraint curves come from Eq. J46t 
and Eq. ( 147b . Although it is not clear from the plot, the limits m — » 0, 
m — * 1 and p — > are all acceptable and yet give no dynamical DE. 
Indeed these are exactly the limits in which we recover standard GR. 



The most interesting part of this parameter space is the limit 
to — > 0, which is one of the limits in which we should recover 
general relativity. The theory then becomes 



f(R) w R - (m 2 + 2A) + p 2 m In (R/p 2 ) 
with the Chameleon-like (singular at <f> = 0) potential 



(55) 



s/6 Mpl 



2A- 



mp' 



In 1 



(56) 

In this limit we are forced to fix (p 2 /2 + A)M p ! 1 = px(0). for 
different values of p. The DE energy equation of state param- 



eter is ojx = — 1 — 0.05mp 2 / Hq. The tightest solar system 
,2 in this limit is from 



constraint on /i z in this limit is from \ j — 1| in Eq. d46*b which 
gives rap? < 6 x 10~ 6 Hq. The equation of state parameter 
for DE is then constrained to be |cjx + 1| 0.3 x 10 -6 which 
is definitely unobservable. 

Finally we note that the ultimate fate of the f(R) 
chameleon is different from that of the original model. This 
is because V(cf>) actually does have a minimum relevant for 
cosmological energy densities. This is due to the <f> depen- 
dence of the AMp\x~ 2 term in Eq. d35l l. which is absent in 
the original models. Eventually will settle into this mini- 
mum and the universe will enter an inflating de Sitter phase, 
much like the fate of a universe with a simple cosmological 
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constant. The original model on the other hand eventually en- 
ters a quintessence like expansion. However, this distinction 
is unobservable today. 

In conclusion to this section, we have found a previously 
unstudied class of f(R) theories that gives acceptable local 
gravity by exploiting the Chameleon effect. For the allowed 
parameters of this model, there is no interesting late-time cos- 
mological behavior (observably dynamic DE). That is not to 
say that these models have no interesting physics — there may 
indeed be some interesting effects of such models for future 
solar system tests ll38ll or on large scale structure ll50ll . and 
this might warrant further study in the context of f(R) mod- 
els. We also noted that the f(R) model is subtly different 
from the original Chameleon model. It does not violate the 
weak equivalence principle , so solar system constraints are 
less stringent and the ultimate fate of the universe is now sim- 
ply an inflating de Sitter spacetime. 

This mechanism might also be a starting point for con- 
structing working modified gravity models which do give non- 
vanilla DE, somehow exploiting this mechanism more effec- 
tively and to bridge the gap in Figure [2] between solar system 
constraints and non-vanilla DE. We suspect they will not be as 
simple as the one presented. This mechanism may also be rel- 
evant for attempting to understand the Newtonian limit of the 
artificially constructed f(R) models that reproduce an exact 
expansion history. We make this claim because an important 
property of the model presented in [28] is that the parameter 
B oc f"(R) is a rapidly growing function of the scale factor 
a. For small f"(R), one can show that the mass curvature 
of V is to| ~ 1/ f"(R). Hence, in this theory the mass of 
the scalar field during cosmological evolution is large at early 
times and small at late times, as in the Chameleon models. 
A more detailed analysis, beyond the scope of this paper, is 
required to see whether non-linear effects play a part in the 
Newtonian limit of these theories. 



IV. MASSIVE f(R) THEORIES 

We now consider arguably more realistic f(R) theories, 
namely polynomials f(R) = -2 A + R + aR 2 + bR 3 .... 
These theories have been extensively studied, especially for 
quadratic f{R)\ see l33ll and references therein. They are 
more natural from the point of view of renormalization and 
effective field theories: a high energy completion of grav- 
ity would allow us to find these higher order terms. How- 
ever, common wisdom would have the higher order terms 
suppressed by inverse powers of M p j and would force us 
to include other terms of the same mass dimension such as 
R^ V R^ V . Despite this, we wish to explore the phenomenol- 
ogy of such polynomial f(R) theories and hence constrain 
them with cosmological observations. In doing so, we will 
explore the full range of values for the coefficients (a, b, ...) of 
the higher order terms to be conservative rather than assume 
that they are order unity in Planck units. 

This class of theories can only match the the currently ob- 
served cosmic acceleration via an explicit cosmological con- 
stant term /(0), giving the identification A = 3H 2 £Ia, so 



there is no hope of dynamical DE. Rather, these theories are 
more relevant to very early universe cosmology where R is 
large, and hence some of our results will be quite speculative. 
Consider for simplicity the two-parameter model 



f(R) = R + R 



XR 



(57) 



We restrict to the parameter range /i 2 > and < A < 1/3, 
so that the resulting potential V has a stable quadratic mini- 
mum and is defined for all <j>. The Einstein frame potential for 
(f) or x is given by 



where 



Ve(x) 



ffS 3A 



2x 2 



q 2 (1 + 2Xq) , 



y/i - 3A(1 - x) - 1 



(58) 



(59) 



is the larger of the two roots of 1 — x + 2q + 3Aq 2 (this ensures 
that the resulting potential has a stable minimum). We plot 
this potential for various A in Figure [3] 



0.2 r 




0/Mpl 



FIG. 3: Potential for the f(R) model in Eq. (57JI with various values 
of A. Notice how the A = case has an asymptotically flat potential 

as <fi — > co. 

We will first explore the possibility that <f> is the inflaton, 
then discuss other constraints from our knowledge about the 
early universe. Figure @] summarizes our constraints. 



A. f(R) inflation 

The possibility that higher order corrections to the gravi- 
tational Lagrangian might be responsible for a de Sitter in- 
flationary period was pointed out early on in the inflationary 
game [ 1 2y . For A = 0, the potential V(<f>) is very flat for 
large <f>, which is perfect for inflation. This model was exten- 
sively studied in |[5lll52ll53ll . which did find a viable inflation- 
ary model. We now search for possible inflationary scenarios 
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of e-foldings of inflation for a mode k is 



lxlO" 5 



FIG. 4: Constraints on the cubic f(R) model. The thin blue/grey 
sliver correspond to observationally allowed f(R) inflationary sce- 
narios. Shaded are regions we may rule out given a measurement of 
the tensor to scalar ratio r and the assumption that they were gen- 
erated by a period of slow roll inflation in the early universe. The 
r — 0.05 and r — 0.01 are the most realistic curve, in the sense that 
future experiments are sensitive to such values as low as r = 0.01 



with A 7^ that are consistent with current observations. This 
question was already considered in [54], which found A <C 1; 
however, we wish to be more quantitative in light of the latest 
CMB measurements. 

As usual in these models, it is important to keep careful 
track of whether we are working in the EF or the JF: recall 
that the potential V is defined in the EF, while matter is most 
naturally considered in the JF. Nonetheless, we will argue that 
the inflationary predictions are exactly the same as those of 
general relativity plus a normal slow rolling scalar field with 
potential V (</>). The argument goes as follows. Slow roll in- 
flation works normally in the EF where the graviton and scalar 
field have canonical actions. In particular, the EF is where one 
should calculate the spectrum of tensor and scalar mode fluc- 
tuations. Re-heating and the transformation of fluctuations 
in 4> to adiabatic density fluctuations also works as usual in 
the EF, because at this time the cosmic fluid is relativistic 
and hence governed by the same equations of motion in both 
frames. After reheating, <fi is frozen out at the minimum of 
V, with (f> = and \ = 1> so there is no longer any distinc- 
tion between the JF and the EF (g^ — g^ u ). Calculations for 
A = were performed both as above and in the JF in 15 311 . 
and the results were found to be consistent as expected. 

Using this idea, calculating the inflationary predictions is 
straightforward. Using Eq. (fTTT i, we can estimate the reheat- 
ing temperature as T RH « 1.3 x 10- 2 g~ 1/4 (N s fi 3 /M pl ) 1/2 , 
where N s is the number of minimally coupled scalar fields 
into which <fr decays (it decays most strongly into these fields). 
Then the scale factor (normalized to a = 1 today) is 



a C nd = 7.5 x 10 



.52 ( M 

Afpi 



-1/6 



-1/12^1/6 



(60) 



at the end of inflation. Integrating the slow roll equations of 
motion, </>' = —V'(<j))/3H, and assuming A <C 1, the number 



3 arctanh(\/A'?) 
2VX 



~ln(l + 2g fc )+iVo(A). (61) 



Here No is a small number defined such that Nk(q en d) = 
at the end of inflation, where q = q enc [ w 1/V3, and q k is 
related to the conformal factor Xk = 1 + 2<?fc + 3Aqf when the 
mode k crosses the horizon: 



H k » ii/y/u = ke N "/a ewi 



(62) 



This particular mode will have a scalar fluctuation amplitude 
(also referred to as 5 2 H in the literature) 



Qk - -,nnn_2. I TJT 
Pi 



12007r 2 e fc I M 



(63) 



where the slow roll parameters (using the definitions in lf55tl ) 
are 



(i -Ml? 



2(1 + Ag|) 
3qk 



(64) 



We then use these to find the the scalar spectral index n s = 
1 — 6e + 2rj, the ratio of tensor to scalar modes r = 16e 
etc. Using the combined WMAP+SDSS measurements [2] 
Q = 1.945 ± 0.05 x 10~ 5 for modes k = 0.002/Mpc we can 
use Eq. doTll63b together to fix ji. For A — * the result is 



(65) 
(66) 
(67) 



H w (3.2 ±0.1) x 10^ 5 M pl , 
n s w 0.964, 
r k, 0.0036, 



which is consistent with the both the theoretical results of 
15311 and recent observational constraints [1,2]. 

In addition, n s is sensitive to the value of A. The observa 
tional constraint 0.937 < n s < 0.969 (68% C.L.) from [20 
translates into a strong upper bound on A: 



A < 4.7 x 10" 



(68) 



This is an example of the usual fine tuning that is needed for 
observationally allowed inflationary potentials and is consis- 
tent with the findings of [54]. More precisely the values of 
H, A appropriate for inflation are shown in Fig|4j 



B. Other constraints 

Above we explored the possibility that <f> was the inflaton. 
Let us now turn to the alternative possibility that <f> is not the 
inflaton, and compute miscellaneous constraints on the param- 
eters n and A when they are varied freely. We will first con- 
sider the fifth force mediated by <j), then investigate how the 
scalar field behaves dynamically in the early universe, where 
the most interesting effect comes from considering a period of 
slow-roll inflation driven by some other scalar field. As noted 
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in Section [TTJ the dynamics of <\> is still governed by an effec- 
tive potential d27l ) which is important when there is a compo- 
nent of matter whose energy-momentum tensor has nonzero 
trace. To begin with, we ignore any effect that such a term 
may have on the minimum of V for these polynomial models, 
which is a good approximation if \Tfi\ <C p 2 M 2 y We will see 
that for the first few constraints that we derive, this will in- 
deed be the case. Then we will return to the question of where 
this is a bad approximation, which will naturally lead to our 
discussion of slow-roll inflation. 



1. Fifth force constraints 

The minimum of the effective potential lies at x = 1, 
(j) = 0. The curvature of this minimum is m| = /i 2 /6. Hence 
we can get around solar system constraints simply by making 
p large enough so that the range of the fifth force will be small. 
Clearly it must have an range smaller than the solar system, 
otherwise, as was discussed above, it will violate the bound 
on the PPN parameter 7. (Recall that there is no Chameleon 
effect here, so A = 1 in Eq. (|43T > and 7 = 1/2.) For smaller 
scales, we consider searches for a fifth force via deviations 
from the inverse square law. The profile for a quadratic poten- 
tial, i.e, Eq. d32l . gives a Yukawa potential between two tests 
masses mi and m-i- 

, . mi m-r e~ m * r 

V ^ = ~ a ^Ml x — (69) 

where a = 1/3. For this a-value, a fifth force is ruled out 
for any Compton wavelength m^ 1 ranging from solar system 
scales down to 0.2mm, where the lower bound comes from 
the Eot-Wash experiments [39]. This bound translates to 

p > 1.0 x l(T 3 eV. (70) 

This implies V(<p) ~ p 2 M 2 y 3* /War, a typical solar system 
density, so for this constraint we were justified in ignoring any 
effects of the density-dependent term on the minimum of V c « ■ 



2. Nucleosynthesis constraints 

Niven this preliminary constraint from local gravity tests, 
let us now consider the cosmology of in the EE We may 
approximate the potential around the minimum by a quadratic 
potential Kff (4>) ~ {p 2 1^)4> 2 , which is valid for \<f>\ < M pl . 
The interesting behavior will come during the radiation domi- 
nated epoch, so in Eq. d24b we take p(a) ss pn(a) oc a~ 4 , and 
we ignore the term in Eq. d25l ) to find the cosmological 
equations of motion 

3H 2 M 2 l = PB{a) + ^+ l -{<p')\ (71) 
4>" + 3H<t>' + ?-cf> = 0, (72) 



where the primes denote d/dt. There are two interesting lim- 
iting behaviors, corresponding to H ^> p and H <C p, which 
we will now explore in turn. 

For H ^> p, the friction term in Eq. (l72l dominates, and <f> 
is frozen out at some value 0* with d<j>/dt = 0. The energy 
density of <p is subdominant in Eq. (TTTT ). Therefore, in the 
EF we have the usual radiation dominated expansion, and in 
the JF using Eq. ( 1231 and Eq. d2Tb we have the same FRW 
expansion with a different effective Newton's constant G* N : 
3H 2 = 8nG* N p(a) oc a~ 4 , where 

For H -C /it, on the other hand, assuming < M p i, the 
field <p starts to oscillate with frequency p/y6 and an am- 
plitude that redshifts as a~ 3 / 2 . Hence in the EF, the energy 
density of cf) in Eq. ( T7TT > from these zero momentum field os- 
cillations is p^ — (p 2 /12)(j) 2 + 4>' 2 /2 w p^(a/a„)^ 3 , where 

P*6 ~ (m 2 /12)0 2 . Mapping back into the JF, and averaging 
over a cycle of this oscillation, we obtain the Friedmann equa- 
tion 

3fT a A^ = p il (o) + |p;(a/o*)- 3 l (74) 

where the additional factor of 1 /2 comes from the averaging 
of the oscillations of G* N in Eq. d73l ), as is discussed in more 
depth in 

The crossover between these two behaviors occurs when H 
is comparable to fi, and given the laboratory tests of gravity 
above we can say that this must occur when the universe has at 
least the temperature T* > ITeV. We were therefore justified 
in assuming radiation domination in our calculation. 

Let us examine further the zero momentum oscillations of 
4> that give this extra non-relativistic energy density. In the 
absence of some mechanism (such as an extra period of low 
scale inflation 15711 ). we expect the initial amplitude of oscil- 
lations to be of the order M p \. This is because the potential 
in Figure [5] varies on the scale of M p \ independently of the 
height of V. Hence in the absence of any other scale, the ini- 
tial amplitude must be around this size. Recall that at the onset 
of oscillations, H ~ /i, so the initial energy density of these 
oscillations is 

p% ~ M p V ~ H 2 M 2 , ~ p R {a*). (75) 

This energy density subsequently grows relative to the radi- 
ation density component, quickly forcing the universe into a 
matter dominated period of expansion. This is unacceptable if 
this component does not decay before the onset of BBN, be- 
cause then at the time of BBN the expansion would be much 
faster than the normal radiation dominated expansion, which 
would be inconsistent with observed primordial abundances 
& 

The fact that <f> interacts weakly with other particles (the 
vertices in Eq. ( fT7b are suppressed by 1 /M p \) so that </> de- 
cays too slowly is exactly what is known as the cosmological 
moduli problem. To be more precise, we can use Eq. ( fTTI i to 
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estimate the decay rate of zero momentum modes into other 
massive particles: 



effects are extremely model dependent, and it is hard to say 
anything definitive at this point. 



E 



96ttM^ 



MJ3847T 



E 



(76) 



where the sums are over minimally coupled scalar particles 
and fermions with masses 2m s , 2m f < m^. The require- 
ment r^, > Hbbn translates into the constraint p > 100 TcV 
for the Standard Model. One would expect the bound on p 
to be slightly smaller if one includes other particles that have 
not been detected yet with mass smaller than lOOTcV. This 
constraint should not be taken too seriously, however, because 
the moduli problem may hypothetically be resolved by elec- 
troweak scale inflation 115811 or even by a brief second period 
of inflation at the electroweak scale 1571. 



3. Density dependent forces 

We now consider how the extra density dependent term in 
V e s may effect cosmology. In other words, when can we not 
neglect the forcing term Tjf of Eq. ( l25b ? After (f> enters the 
oscillating phase when /i> H, the extra term has little effect 
on the minimum since then it is small compared to the size 
of the potential itself (V ~ p 2 M^). As a result, <p simply 
oscillates as expected. Before the crossover, when <j> is frozen, 
we showed that the universe must be radiation dominated so 
that in particular, as T£ <C p = SH 2 M 2 X during this phase, 
the Hubble friction will dominate compared to the force term 
of Tjf in Eq. d25l ), and we were justified in claiming that cf) 
is frozen out. The cosmology here does not suffer from the 
instability that plagues Eq. ([T). 

There are, however, some exceptions that might lead to in- 
teresting constraints. First, consider a relativistic component 
i of the cosmological plasma that becomes non-relativistic 
and dumps its energy into the other relativistic components. 
In this case, — T£ ~ (9i/g*)p f° r a period of about one e- 
folding, so <j> receives a kick and is displaced by an amount 
A<p ss (.9i/<7*)Mpi/v / 6 l47ll . This might lead to an interest- 
ing effect such as <f> being kicked out of the basin of attraction 
of V. The extreme case would be that (f> does not end up os- 
cillating around the minimum as expected when H ~ p, but 
instead rolling down the tail of V, an effect which is clearly 
only possible for A ^ 0. In principle, such kicks could even 
invalidate the predictions of BBN: near the onset of BBN, 
annihilation occurs, displacing <fi and consequently changing 
significantly as per Eq. (l73l l. However, we have already 
shown that tj> must be in the oscillatory phase long before the 
onset of BBN, and we have argued that these kicks have no 
effect while is in the oscillator phase, so in fact this effect 
is unlikely to have relevance for BBN. Such kicks may ef- 
fect other important cosmological dynamics at temperatures 
higher than T > ITeV, such as baryogenesis. However, the 



4. Non-inflation 

Another situation when we cannot ignore the density de- 
pendent force on <fr is during inflation. Here is large for 
many e-foldings. Remember that in this section, we are not 
considering </> as our inflaton; instead we consider a slow roll 
inflationary period driven by some other scalar field ip defined 
in the JF. We wish to examine the effect a modified gravity 
Lagrangian such as Eq. (l57l i has on the inflationary scenario. 
In particular, we will be interested in situations where infla- 
tion by the field ip does not work, being effectively sabotaged 
by (p. We will discuss the generality of these assumptions at 
the end. 

Such models have been considered before in the context of 
both the A = models II591 16CH. 16111 and other generalized 
gravity models ll62ll . There the goal was generally to make the 
inflationary predictions more successful, focusing on working 
models. 

In the JF, consider a scalar field ip with a potential U(ip). 
We assume that ip is slow rolling; dip/dt S3 —U'(ip)/3H(t). 
This is the assumption that 



dhp_ 
dt 2 



~dt 



(77) 



which must be checked for self-consistency once we have 
solved for H(t). We can now easily calculate H(t) by first 
working in the EF and mapping back to the JF. The equations 
of motion in the EF Eq. (1241125b become 



3iJ 2 M 2 = -cf,' 2 + V eS (<f>) , <P" + 3H4? = -Kff W (78) 



Kff (0) = V(cp) + t/(V0x~ 



(79) 



It is interesting that a constant vacuum term in the JF does not 
translate into to a constant term in the EF. See Figure [5] for 
some examples of the effective potential V e s; we see that for 
large enough U(tp) ^S> p 2 M 2 Y , the minimum vanishes. One 
finds that there is no minimum of the effective potential for 
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(80) 



In particular, there is always a minimum for A = 0. 

The resulting behavior of the inflaton i/j depends on the size 
of U(tp) compared to p 2 M 2 v For small U(ip) <§C p 2 M 2 x , it is 
clear that there is a stable minimum around which (p will os- 
cillate. In this situation, the effective potential has a minimum 
at cf) w with value approximately V e s{4> ~ 0) w U(ip), 
so after the energy density of 4> oscillations redshift away, we 
are left with an exponentially expanding universe with x ~ 1> 
S&M^ w U(tp) and H S3 H. Hence in the JF, gravity 
behaves as it normally would in general relativity: for a flat 
potential, the slow roll conditions are satisfied, and inflation 
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2 4 6 8 

FIG. 5: Effective potential for the polynomial model Eq. d57t . with 
various JF inflationary energy densities u(ip), here A = 0.1. 



driven by ip works as it normally would. This is the expected 
situation, and will happen for fi ps M p \. 

On the other hand, we now show that when U(ip) 3> [^M 2 ^ 
and when there is no minimum of the effective potential (A ^ 
0), we get a contradiction to the assumption that ip was slow 
rolling. Hence we show that it is not possible for ip to drive 
slow-roll inflation. For large U(ip) 3> fi 2 M^, the potential 
may be approximated as V e g p=> U(ip)x~ 2 - We treat U(ip) as 
a constant and find that there is an exact attractor solution to 
Eq. d78l ) of the form \ ~ t and 5 ~ i 3 / 4 . Mapping this into the 
EF, we find the behavior a ~ i 1 / 2 , i.e., a period of radiation 
dominated expansion analogous to the 0MDE of [23] More 
specifically, we find 

SM^H 2 re U(ip)a- 4 (81) 

This is clearly not an inflating universe. So the slow roll as- 
sumptions of Eq. ( TTTb are not consistent in this case. We there- 
fore conclude that it is not possible for ip to drive slow-roll 
inflation. 

Instead, ip dumps most of its energy U(ipo) into radia- 
tion, and as before, <p is left frozen at some point 0* until 
t/(i/ , o)a -4 ~ (J?M^. After this, <p can either drive an in- 
flationary period itself as in the original discussion of f(R) 
inflation, or if 0* is not in the basin of attraction, it will roll 
down the tail of V. In neither situation has ip inflated our uni- 
verse. From this combination of inflaton ip plus f(R) gravity, 
we only get satisfactory inflation if (fj,, A) lies in the region of 
parameter space appropriate for f(R) inflation (the blue/grey 
sliver in Figure|Ui and if <p* sits at a point which allows for the 
required number of e-foldings. 



5. Gravitational wave constraints 

It is well-known that inflation produces horizon-scale grav- 
itational waves of amplitude Q t ~ H/M p i, so that the energy 
scale of inflation can be bounded from above by the current 



observational upper limit Q t < 0.6 x 1CT 5 d H and per- 
haps measured by a detection of the gravitational wave signal 
with future CMB experiments [49]. However, because of the 
EF-JF duality, one needs to carefully define what one means 
by "the energy scale of inflation". The bound from the above 
argument simply precludes inflatons with a given energy den- 
sity U (ip) in the JF, but U(ip) is merely a parameter which 
does not necessarily set the energy scale of inflation. In ad- 
dition, we cannot derive a constraint for A = 0, because then 
there is always a minimum in the effective potential and it is 
always possible for ip to slow roll (this situation is described 
in greater depth in II591I60I1 ). 

To make these ideas more concrete and resolve both of 
these ambiguities, we will operationally define the energy 
scale of inflation to be the one that makes the standard GR 
formula for the gravitational wave amplitude valid. It is clear 
that the amplitude of gravity waves should be calculated in the 
EF where the metric has a canonical action. The result is then 
passed trivially into the JF after inflation and when (p = 0. 
The Hubble scale H then sets the size of the fluctuations, but 
it is a complicated model dependent calculation to find exactly 
when the relevant fluctuations are generated. However, there 
is a limit to the size of H for which the EF is approximately 
inflating, and so gravity waves are being generated. Follow- 
ing the discussion above of non-working inflatons, we demand 
that <p must be slow rolling down the effective potential V c g 
defined in Eq. d79l ), for both frames to be inflating. In this 
situation, both scalar and gravity modes are being generated. 

The procedure is thus to find the maximum value of H (that 
is, the maximum value of V e s) such that <p is slow rolling. We 
then maximize this H with respect to the parameter U(ip) to 
find the largest amplitude of gravity waves that can possibly 
be produced. At each step in this procedure, we wish to be 
as conservative as possible; for example, we define slow roll 
through the slow roll parameter constraints e < 1 and \rj\ < 2 
to allows for the possibility of power law inflation. 

As an example, consider the A = case. Here it is possible 
to show that for <p to be slow rolling, it must satisfy 

where </> sr always lies to the left of the minimum of V c s. The 
maximum Hubble scale in the EF for a given U(ip) is then 

H 2 < max{V r cff (0 sr )/3Af 2 1 , // 2 /24}. This is maximized 

for large U(ip)/ fJ?M 2 b with the result that H 2 < ^ 2 /6. 
This translates into a constraint on the maximum gravitational 
wave amplitude 

Qf AX « 0.04-^ r MAX « 5 x lO 6 -^- (83) 

* Mpi M 2 

Given a measurement of the tensor to scalar ratio r, this places 
a limit on /i: 

\l > 3 x KTV /2 M pl (84) 

Numerically, we find similar results for non-zero A. We plot 
examples of this constraint in Figure |4] combined with the 
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already discussed working f(R) inflationary models. Note 
that for a given r, it is important that this constraint lies below 
the corresponding working inflationary model with the same 
r; fortunately, it does. 

If gravitational waves are not detected, then this argument 
gives no lower bound on /i. In particular, it is possible that 
inflation occurred at the electroweak scale, in which case the 
constraint // > 2 x 10 _3 eV is the best we can do. 

Note that we completely ignore the production of scalar 
fluctuation modes for this argument. This is because the scalar 
modes are much more difficult to calculate, since there are two 
scalar fields in the mix, tp and <fi, which are canonically de- 
fined in different frames. But the scalar modes are also model 
dependent and one should generally be able to fine tune U (ip) 
to give the correct amplitude and spectral index without af- 
fecting the above argument. This more complicated problem 
was considered for chaotic inflation with R 2 gravity in [61]. 

This constraint applies only to slow-roll inflation models. 
There are classes of fast-roll inflation, but these models have 
problems of their own and generally fail to reproduce the re- 
quired scale invariance (see B63I1 for a review). 

Finally, let us discuss some inflaton models that might cir- 
cumvent this constraint. It is possible to add an inflaton in the 
EE However, this theory is then not conformally equivalent to 
an f(R) theory: the two scalar fields tp and <p get mixed up. 
Hence it is not in the class of models we set out to constrain. 
Another possibility is to add an inflaton which is conformally 
coupled to gravity and with a V oc ip 4 potential. This does 
not change from frame to frame and so inflation might be 
expected to work. However, it was shown by fl64ll that non- 
minimally coupled scalar fields cannot drive inflation. 

In any case, if gravitational waves are found, then this con- 
straint must be thought about seriously when using such f(R) 
models in other astrophysical or local gravity situations. 

V. CONCLUSIONS 

We have have searched for viable f(R) theories using the 
wealth of knowledge on scalar scalar tensor theories to which 
f(R) theories are equivalent. We studied two classes of mod- 
els: the f(R) Chameleon and massive f(R) theories, which 
may well be the only classes of models that can be made con- 
sistent with local gravity observations. 

The f(R) Chameleon that was studied is a special kind of 
scalar field which hides itself from solar system tests of grav- 
ity using non-linear effects associated with the all-important 



density dependent effective potential. It was shown that de- 
spite this Chameleon behavior, solar system tests still preclude 
the possibility of observably dynamical DE; the best we could 
do was \wx — 1| ^ 0.3 x 10~ 6 for the effective DE equation 
of state parameter wx relevant for the dynamics of the expan- 
sion. There are of course interesting effects of the Chameleon 
both for local gravity ll37ll and on cosmological density per- 
turbations [50], and these may be worth future studies in the 
context of f(R) theories. 

The massive theories were found to be more relevant for 
very high energy cosmology, so the conclusions were more 
speculative. First, the scalar field may be the inflaton, in which 
case we found the required polynomial f(R) to be quite fine 
tuned as is usual for inflationary potentials. If the scalar field 
was not the inflaton, then we saw that possible instabilities 
could spoil both inflation and Big Bang nucleosynthesis, giv- 
ing interesting constraints on the shape of f(R). If primor- 
dial gravitational waves are detected using the CMB, then the 
most naive models of inflation have serious problems unless 
the mass of the /(i?)-scalar is very large; a measured scalar to 
tensor ratio of r = 0.05 requires /j>7x 10~ 5 M p i. If gravity 
waves are not found, then the best we can say comes from the 
Eot-Wash laboratory experiments constraining the extent of a 
5th force; ^ > 2 x lCT 3 eV. 

General relativity adorned with nothing but a cosmological 
constant, i.e., f(R) = R—2A, is a remarkable successful the- 
ory. As we have discussed, a host of observational data prob- 
ing scales from 10~ 2 m to 10 26 m not only agree beautifully 
with GR, but also place sharp constraints on the parametrized 
departures from GR that we have explored. In particular, both 
viable classes of f(R) theories that we studied were found 
to have no relevance for dynamic dark energy that is observa- 
tionally distinguishable from "vanilla" dark energy, i.e., a cos- 
mological constant. Since we have no good reason to believe 
that there are additional viable classes of f(R) -theories, it ap- 
pears likely that no viable f(R) theories can produce the sort 
of interesting non-vanilla dark energy that many observers are 
hoping to find. 
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